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ABSTRACT 

We derive the stacked 1.4 GHz flux from FIRST (Faint Images of the Radio Sky at Twenty Centimeters) 
survey for 811 K+A galaxies selected from the SDSS DR7. For these objects we find a mean flux density of 
56 ± 9 /iJy. A similar stack of radio-quiet white dwarfs yields an upper limit of 43 /iJy at a 5(T significance to the 
flux in blank regions of the sky. This implies an average star formation rate of 1.6 ± 0.3 year~^ for K+A 
galaxies. However the majority of the signal comes from ~4% of K+A fields that have aperture fluxes above 
the 5(7 noise level of the FIRST survey. A stack of the remaining galaxies shows little residual flux consistent 
with an upper limit on star formation of 1.3 Mq year~^. Even for a subset of 456 'young' (spectral ages < 250 
Myr) K+A galaxies we find that the stacked 1.4 GHz flux is consistent with no current star formation. Our 
data suggest that the original starburst has been terminated in the majority of K+A galaxies, but that this may 
represent part of a duty cycle where a fraction of these galaxies may be active at a given moment with dusty 
starbursts and AGNs being present. 

Subject headings: galaxies: starburst — galaxies: evolution — radio continuum: galaxies 



1. INTRODUCTION 

K+A galaxies, also known as post-starburst galaxies (PSG), 
show spectra of a strong B aimer absorption series super- 
posed over a K-giant dominated spectrum typical of early- 
type galaxies, implying the recent termination (< 1 Gyr) of a 
significant episode of star formation in an otherwise quiescent 
stella r population (Dressier & Gunnll 19831: ICouch & SharplesI 
119871) . These objects may be the best examples of galax- 
ies transitioning from the 'blue cloud' (of star-forming ob- 
jects) to the 'red sequence', and have often been iden- 
tified as p ossible progenitor s of t he lenticular population 
in clusters (j Yang et al. 2004 120081) , rejuvenated early-type 
galaxies (e.g., Panuzz o et al.l 1200% or the descendants of 
the blue galaxies observed in intermediate redshift clusters 
(IZabludoffetal.llT996l: iPoggianti et an[l999h . 

The mechanism by which star formation is initiated and/or 
quenched in these galaxies is still unclear. The majority of 
local K+A galaxies lie in the general field rather than clus- 
ters ( Blake etal.ll2QQl iGotol 1200 71: I Vergani et al.l IMoh and 
the processes that trigger and halt star formation may differ in 
these environments. At the same time, in the intermediate red- 
shift clusters where K+A galaxies were originally discovered, 
a significant f raction of current star for mation is obscured 
by dust (Due et al.ll2002l: lSaintonge et"an[2008: Dressier et al.l 
[2009; Haines et al. 200^. Dusty starbursts tend to have strong 
B aimer absorption but only weak Oil emission as the young 
stars and HII regions that produce the emission lines tend to 
lie in regions of high extinction while the longer lived A- 
stars can migrate out of their native molecular cloud: this 
would be a viabl e model to explain the K+A spectral class 
(iPoggianti & Wul 12000) . Therefore it is still unclear whether 
K+A galaxies are truly "red and de ad," although the 24/im 
observations of "Dressie r et al.l (120091) in Abell 851 imply that 
star formation has largely ceased in these objects. 
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Obscured star formation may also be detected via radio 
emission at 20 cm (1.4 GHz), which is produced by syn- 
chrotron radiation from high energy cosmic rays originating 
in supernova sh ells and y ields a measure of the massive star 
formation rate (iCondonI 119921: iKennicuttI 119981) . Although 
iSmail et aH (1199 9^ found evidence of recent star formation 
in 5 K+A galaxies, other studies find little evidence of on- 
going star formation within small samples ("Mi ller & OwenI 
120011: iGoto 2004). However, Buvle et al. (200j) found that 
most K+A have substantial gas reservoirs, similar to those of 
spirals of the same luminosity, and suggest that K+A galaxies 
may be observed during a hiatus in an episodic star formation 
history or that current star formation may be obscured. 

In this Letter we use a stack of radio images to further test 
the current star formation rate and/or AGN activity in spec- 
troscopically selected K+A galaxies. A description of the 
data and the analysis is provided in the next section, while 
we interpret and discuss our results in section 3. We adopt 
the latest WMAP cosmological parameters with Q^m = 0.27, 
VtA = 0.73 and Hq = 71 km s-^ Mpc-^ 

2. FIRST RADIO STACKING OF THE SDSS DRV K-fA SAMPLE 

We select a sample of 811 K+A galaxies f rom the Data Re- 
lease 7 of the SDSS (lYork etaLjllO OO: Aba zaiian et al.ll2009h 
using an updated catalog of'Goto (2007). Only objects clas- 
sified as galaxies with a spectroscopic signal-to-noise > 10 
per pixel are considered. The selection criteria of K+A galax- 
ies are equivalent widths of Ha > — 3.0A, H(5 > 5.0A and 
[Oil] > -2.5 A, where emission lines are negative. Galaxies 
of redshift 0.35 < z < 0.37 are excluded from the sample 
due to the 5577A sky feature. The selected sample of K+A 
galaxies have redshifts ranging 0.02 < z < 0.4. 

We then use data from the FIRST survey (iBecker et al.l 
IT995I) to derive a mean radio image of K+A galaxies. Typ- 
ical detection limits for a single FIRST image are ~1 mJy, 
which, for the mean redshift from our SDSS DR7 sample of 
K+A galaxies, would only allow us to detect star formation 
rates in excess of 30 Mq yr~^. 

For each of the 811 K+A galaxies covered by the survey 
we cut out a V square from the FIRST database. We use the 
median stacking method of .White et al.. (2007.) to create the 
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resulting image of 8 1 1 K+ A galaxies shown in the left panel 
of Fig.[Tl We then measure the flux from the combined image 
in an aperture equivalent to three FIRST beams from which 
we derive an average K+A flux of 56 ± 9 /iJy. 

To assess the significance of this result, we need to estimate 
the level of noise in the image. As FIRST is a survey, the 
image cutouts may not be fully cleaned of artifacts, such as 
side-lobes from distant radio sources. To do this, we use a 
Monte Carlo simulation with a sample of known radio-quiet 
objects, 8495 white dwarfs from the SDSS (Eisenstein et al. 
^006). We create 10,000 stacks of 811 randomly selected 
white dwarfs (equivalent to the stack of K+A galaxies) and 
measure the flux in the same fashion to derive the mean flux 
from a supposedly radio-quiet sample. 

This procedure allows us to estimate the level of noise 
present in the stacked image and therefore to determine the 
significance of our detection for K+A galaxies. A histogram 
of the fluxes from our Monte Carlo simulation is shown in 
Fig. [21 Since we know that these sources are radio quiet, we 
can estimate a 5 a significance flux for sources to be consid- 
ered real. We find this flux to be 43 /iJy. Our K+A stack is 
found to have a significant detection above the 5 a noise level. 

We can now compute the absolute luminosity at 1.4 GHz 
using our chosen cosmology and the average redshift of the 
K+A galaxies (z = 0.14). 

LlAGHz = ^TTDlSiAGHzil + ^)"/(l + z), (1) 

where Dl is the luminosity distance, Siaghz is the flux den- 
sity, (1 + z)" is t he color correction an d 1/(1 + z) the band- 
width correction (Morrison et al.ll2003 l). We assume that the 
radio emission is dominated by synchrotron radiation such 
that S (X (Condon 1992). This yields our measurement 

of radio power, which we can convert to a SFR using 

SFR (Mq yr-i) = 5.9 x IQ-^^LiaghzC^ Hz-^), (2) 

from lYun et al.l (|2001|) . which assumes a Salpeter initial mass 
function between 0.1 and 100 solar masses. This yields a star 
formation rate of 1.6 ± 0.3 M© yr~^ for our average K+A 
galaxy at < z >= 0.14. 

However 79 of our sources are measured to have 1.4 GHz 
fluxes above the 3 a noise level of the FIRST survey. The 
FIRST beam has an RMS of 0. 15 mJy and we find 79 galaxies 
(~10%) have aperture fluxes in excess of 3a (450 /iJy) within 
an aperture of 3 beamsizes, 31 of which have fluxes in excess 
of 5cr (750 /iJy). 

In Fig. [3] we show the redshift-luminosity distribution of the 
targets with aperture fluxes above the 3a and 5a limits. Visual 
inspection of these individual outlier frames have shown that 
for the 31 galaxies above the 5a limit, about 80% look like 
clear individual detections (i.e., centrally concentrated, likely 
to be associated with our optical target, not in a frame that 
is significantly noisy), while for the 48 that are above the 3a 
limit but less than the 5a limit, only about 10% of these are 
clear detections. Overall only about 4% of our sample of 81 1 
galaxies show evidence of significant ongoing radio activity, 
either from star formation or AGN activity. 

Removing all sources with measured fluxes in excess of 5a, 
we create a subsample of 780 galaxies. Stacking this subsam- 
ple yields a mean flux of 36 /iJy, which is well below the 5 a 
detection limit of 47 /iJy found from a Monte Carlo simula- 
tion creating stacks of 780 white dwarfs. We can place an 
upper limit on the SFR of 1.3 M© yr~^ for this subsample 
with < z >= 0.14. 



When star formation ceases in a galaxy, we expect it to de- 
crease exponentially to low levels, rather than an abrupt trun- 
cation, unless some 'catastrophic' event has removed all the 
available fuel or ionized the gas and prevented further star 
formation. Evidence for a rapid shutdown of star formation 
in these galaxies has been presented by Brown et al. (20Q9|). 
We therefore carry out the same analysis on a subsample of 
'young' K+A galaxies and attempt to detect residual star for- 
mation in these objects. 

In Fig. H we plot US and D4000 from SDSS spec- 
t ra and ove rplot models calculated with GALAXEV 
(iBruzual & Chariot h 2003). We plot the galaxies with signifi- 
cant 1.4 GHz detections as red squares, the rest as gray dots. 
We adopt a Salpeter initial mass function (IMF) and solar 
metallicity as initial conditions. The model galaxies evolved 
over 10 Gyr with an exponentially decreasing star formation 
rate (r = I Gyr). At 10 Gyr, we added an instantaneous star- 
burst (delta function) of mass 1,5, 10, 30 and 50% (relative to 
the old stellar population) after which the SFR returns to zero. 
The dotted lines demarcate the values of the spectral indices 
observed 30, 50, 100, 250 and 500 Myr after the bur st. These 
are the same models used by iYagi et al.i (|2006)) and lGoto et al.l 
(2008). 

We then select the 456 K+A galaxies with burst ages less 
than 250 Myr and an average redshift of < z >= 0.13 and 
stack these galaxies in the same fashion as the complete sam- 
ple. The stacked image is shown in the right panel of Fig. [T] 
from which we measure a mean flux of 61 /iJy. The Monte 
Carlo simulation with white dwarfs gives a 5 a detection limit 
of 59 /iJy. Our younger sample of K+A galaxies is found to 
be at the ba detection limit, giving a SFR of 1.5 ± 0.3 M© 
yr~^. This is consistent with a rapid decline of st ar formation 
in these galaxies as shown by lBrown et al.l (|2009|) . 

3. DISCUSSION 

We have found an average star formation rate of only 
1.6 ± 0.3 M0 yr~^ from stacking radio observations of 811 
general field K+A galaxies in the local universe. 

However, much of this signal appears to originate from 
^4% of active galaxies. Based on the definition of 
ISadler et al.l (12002), which requires radio power above 10^^ 
W Hz~^ and the absence of emission lines, this sample is ap- 
proximately equally split between star forming galaxies and 
AGN (with the latter more prominent at higher redshifts be- 
cause of selection effects). 

For the remainder of the sample, we find an upper limit on 
the SFR of 1.4 M© yr~^. Even a subsample of spectroscopi- 
cally young galaxies does not show a significant detection of 
residual star formation. 

Ou r resul ts are in good agreement with previous work. 
iGotol (|2004|) did not detect evidence of star formation from 
VLA observation s in 36 galaxies drawn from the SDSS DRl 
(iGoto et al.l2003 ) but he was able to set an upper limit of < 15 
M^yr"^ for 15 of the nearest K+A galaxies. Mil ler & OwenI 
(1200 1|) found that only 2 out of 15 K+A galaxies in his sam- 
ple show signs of obscured star formation. Nevertheless, 
iMelnick & de ProprisI (120 12h find that most K+A galaxies in 
this sample exhibit significant excess above the predicted stel- 
lar component at A > 5 fim from WISE data, suggesting that 
dust heated by some unknown source is present in these ob- 
jects. 

iBuyle et aP (120061) pointed out that K+A galaxies contain 
significant amounts of gas and therefore the current quies cent 
star formation may only be temporary, while Poggianti & Wul 
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Fig. 1 . — Left: The median stack of FIRST images of 811 K+A galaxies selected from the SDSS. The image is V square with grayscale units in /iJy. Right: 
The median stack of 427 K+A galaxies with less than than 250 Myr since the end of the starburst. The image is 1' square with grayscale units in /iJy. 
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Fig. 2. — Histogram of fluxes from Monte Carlo stacks of white dwarfs in 
the SDSS. Each stack is equivalent to the stack of K+A galaxies used in the 
analysis and contains 811 white dwarfs images. We find a 5a noise level of 
43 /iJy. 



(I2QQ0I) proposed a model where star formatio n continues in 
K+A galaxies but is hidden by dust. However. iDressler et alJ 
(f2009) detected no 24 /am emission for K+A galaxies in Abell 
851. 

Another possible contribution to the radio flux is that K+A 
galaxies may host an AG N. Most K+A galaxies have signif- 
icant bulges ( Yang et al.l 1200 8) and ther efore should contain 
a supermassive central black hole (e.g., iFerrarese & Merritfl 
|2000). In the general field, K+A galaxies are often involved 
in mergers and interactions (.Blak e et al. 2004; Goto 2005; 
I Yang et al.ll2008h and sh ow inverted color gradients indicative 
of central star formation (lYang et al."2008'). which is expected 
if mergers d rive gas to the center (Hernquistl ll989l) . 

iLiu et al.l (|2007) find evidence of a weak AGN in a nearby 
K+A galaxy. However, .Shin et al.. C2011.) find only a few 



Fig. 3. — We have plotted K+A galaxies from our sample with aperture 
fluxes in excess of 3cr and 5cr as green squares and red circles, respectively. 
The dotted and dot-dashed lines indicates the 3a and 5a FIRST survey limits, 
respectively, as a function of redshift. 

AGNs in a subset of our sample and they argue that these 
active nuclei are not related to the quenching of the previ- 
ous starburst. Brown et al. (2009) argue that AGN feedback 
is only important among the more massive K+A galaxies. 
Vergani et al. (2010) find no evidence for AGN in their sam- 
ple of COSMOS K+A galaxies using deep VLA and XMM 
imaging of the COSMOS field and only a small number of 
star forming objects, in agreement with our results. On the 
other hand Georgakakis et al.' ("2008) argued for a weak corre- 
la tion between AGN and K+A galaxies in the EGS field. 

ISnyder et^.l (12011 ) present a model where K+A galax- 
ies originate from gas-rich mergers and the duration of the 
B aimer-line strong phase is shorter than the commonly as- 
sumed ~1 Gyr by about a factor of 3. In this framework the 
~4% of 'active' galaxies that we find may be the tail end of 
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Fig. 4. — Plot of ¥16 and D4000 for K+A galaxies from SDSS spectra and models. The gray circles and red squares are the measured values for K+A galaxies. 
The red squares indicate galaxies with significant 1.4 GHz detections. The solid lines are the models of an old 10 Gyr exponentially decaying stellar population 
with a burst of 1% (blue), 5% (purple), 10% (cyan), 30% (green) and 50% (yellow) mass relative to the older stellar population. The dotted lines indicate the 
values of the spectral indices when observed 30, 50, 100, 250 and 500 Myr after the burst. 



the distribution, while other objects are already heading to- 
wards the red sequence. In this case, this might yield a con- 
straint on the activity timescale, assuming a 0.3 Gyr duration 
for the K+A phenomenon. If the 4% of active galaxies rep- 
resent the last gasp of activity before the AGN returns to a 
dormant state, the period of visible QSO activity may be es- 
timated to be about 20 Myr, which is broadly consistent with 
current estimates (e.g., see review by MartinildlOOi)). 

Another plausible explanation is that we are actually ob- 
serving galaxies undergoing a series of starb ursts and feed- 
back e pisodes, as in the cycle postulated by iHopkins et al.l 
(120081) to account for the tightness of the red sequence within 
a hierarchical formation framework. In which case we would 
be observing galaxies in the 'on' and 'off phase of such a 
c ycle. 

iShioya et"an (12002 1) present a model where truncated spirals 
evolve through the e(a) (galaxy with strong US and modest 
on emission), a+k and finally k+a phase before settling on the 
red sequence. If we are observing galaxies at random along 
their spectrophotometric transformation, and given the ~ 0.3 
Gyr duration of the most visible k+a phase, the timescale for 
evolution of these galaxies on to the red sequence is long, ~ 7 
Gyr, which is consistent with the above simulations. Only low 



mass objects can then truly evolve on to the red sequence in 
reasonable times. 

A series of questions remain unanswered for future stud- 
ies. Do K+A galaxies host on-going dusty starbursts? This 
appears unlikely in the light of our res ults, but the detection 
of pow erful far infrared excesses by iMelnick & de Proprlsl 
(120121) points to the existence of extra sources of flux that 
are heavily obscured. What is the role of the AGN (if any) 
in modulating the rapid onset and decline of star formation? 
Do these objects contain cold gas and possibly re-initiate star 
formation? Future papers by our group will attempt to address 
some of these issues. 
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